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Abstract 

In this article, new results are presented for the zero-temperature ground-state properties of the 
spin-half transverse Ising model on various lattices using three different approximate techniques. 
These are, respectively, the coupled cluster method, the correlated basis function method, and 
the variational quantum Monte Carlo method. The methods, at different levels of approximation, 
are used to study the ground-state properties of these systems, and the results are found to be 
in excellent agreement both with each other and with results of exact calculations for the linear 
chain and results of exact cumulant series expansions for lattices of higher spatial dimension. 
The different techniques used are compared and contrasted in the light of these results, and the 
constructions of the approximate ground-state wave functions are especially discussed. 
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I. INTRODUCTION 



Two of the most versatile and most accurate semi-analytical formalisms of microscopic 
quantum many-body theory (QMBT) are the coupled cluster method [jXJ, |2], |3|, |], |, ^ |7], || 
and the correlated basis function (CBF) method [|, [[g, [TT], |I§ 0, [16], |T7|, |T|, 
p~9|1 . In recent years such QMBT methods, together with various quantum Monte Carlo 
(QMC) techniques, have been applied with a great deal of success to lattice quantum spin 
systems at zero temperature. Some typical recent examples of such applications include 
Refs. J2U], £T[ [22], [23J [?§ |2§ |2Bj, for the CCM, Refs. |2g 0, for the CBF 
method, and Refs. [33], [35|, |3_7], 3J| [3J5, for the various QMC techniques. Current 
state of the art is such that these methods are sufficiently accurate to describe the various 
quantum phase transitions between the states of different quantum order that exist in such 
abundance for spin-lattice systems. However, each of the above methods is characterised by 
its own strengths and weaknesses. Hence, a fuller and more complete understanding of such 
strongly interacting systems as the lattice quantum spin systems is expected to be given 
by the application of a range of such techniques than by the single application of any one 
of them. In this article we wish to apply the CCM, the CBF method, and the variational 
quantum Monte Carlo (VQMC) method to the spin-half transverse Ising model (for reviews 



of this model see, for example, Refs. [0], ^5], ^6|). The Hamiltonian for this 

model on a lattice of N sites, each of which has z nearest-neighbours, is given by 



H 



(1) 



where the cr-operators are the usual Pauli spin operators and (i, j) indicates that each of 
the zN/2 nearest-neighbour bonds on the lattice is counted once only. We work in the 
thermodynamic limit where N — > oo. We note that this model has an exact solution in 
one dimension [flfifl; and approximate techniques, such as the random phase approximation 
(RPA) |47|] and exact cumulant series expansions PH|, f29| , have also been applied to it for 
lattices of higher spatial dimensionality. For A > 0, we note furthermore that the model 
contains two distinct phases, with a critical coupling strength A c depending on lattice type 
and dimensionality. For A < A c there is non-zero spin ordering in the ^-direction, and hence 
this regime will be referred to here as the ferromagnetic regime. By contrast, for A > A c , 
all of the ferromagnetic ordering is destroyed, and the classical behaviour of these systems 
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is that the spins lie along the positive x-axis. Hence, the A > A c regime will be referred to 
here as the paramagnetic regime. In Sec. II the technical aspects of applying the CCM, the 
CBF, and the VQMC methods to the spin-half transverse Ising model are presented, and in 
Sec. Ill the results of these calculations are discussed. Finally, the conclusions are given in 
Sec. IV. 

II. QUANTUM MANY-BODY TECHNIQUES 
A. The Coupled Cluster Method (CCM) 

In this section, we firstly describe the general CCM formalism |L], |2], |3|, £|, ||, ||, |7|, |[, and 
then proceed to apply it to the specific case of the spin-half transverse Ising model. The 
exact ket and bra ground-state energy eigenvectors, and of a many-body system 
described by a Hamiltonian H, 

H\V) = E g \*); (V\H = E g (*\, (2) 

are parametrised within the single-reference CCM as follows: 

|*> = e 5 |$> ; S = J2SiCt, 

(§| = ($|£e- 5 ; S=l + Y,Si C I- (3) 

The single model or reference state |$) is required to have the property of being a cyclic 
vector with respect to two well-defined Abelian subalgebras of multi-configurational creation 
operators {Cf} and their Hermitian-adjoint destruction counterparts {Cj = (CfY}. Thus, 
|<3>) plays the role of a vacuum state with respect to a suitable set of (mutually commuting) 
many-body creation operators {C/}, 

C7|$> = ,1/0, (4) 

with Cq = 1, the identity operator. These operators are complete in the many-body Hilbert 
(or Fock) space, 

i = |$)($| + ^c+|$)($|cy . (5) 

Also, the correlation operator S is decomposed entirely in terms of these creation operators 
{Cf}, which, when acting on the model state ({(7/1$)}), create excitations from it. We 
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note that although the manifest Hermiticity, ((&\^ = \^) / (^\^)), is lost, the intermediate 
normalisation condition (^l^) = ( < £| l I / ) = (^|$) = 1 is explicitly imposed. The correlation 
coefficients {Si, Si} are regarded as being independent variables, even though formally we 
have the relation, 

(*\s = J;f' t * . (6) 

($|e 5T e 5 |$) 

The full set {Si, Si} thus provides a complete description of the ground state. For instance, 
an arbitrary operator A will have a ground-state expectation value given as, 

A = = (<S>\Se- s Ae s \$) = A({Si,Si}) . (7) 

We note that the exponentiated form of the ground-state CCM parametrisation of Eq. (|3p 
ensures the correct counting of the independent and excited correlated many-body clusters 
with respect to |<&) which are present in the exact ground state |^}. It also ensures the 
exact incorporation of the Goldstone linked-cluster theorem, which itself guarantees the 
size-extensivity of all relevant extensive physical quantities. 

The determination of the correlation coefficients {57,5/} is achieved by taking appropri- 
ate projections onto the ground-state Schrodinger equations of Eq. (H). Equivalently, they 
may be determined variationally by requiring the ground-state energy expectation functional 
H({Si, Si}), defined as in Eq. (|7|), to be stationary with respect to variations in each of the 
(independent) variables of the full set. We thereby easily derive the following coupled set of 
equations, 

5H/5Si = ($|<7fe- 5 #e 5 |$) = 0, / ^ ; (8) 
5H/5Si = (<5>\Se- s [H,Cf}e s \<5>) = 0, 1^0 . (9) 

Equation ([8]) also shows that the ground-state energy at the stationary point has the simple 
form 

E g = E g ({Si}) = ($|e-^e 5 |$) . (10) 

It is important to realize that this (bi-)variational formulation does not lead to an upper 
bound for E g when the summations for S and S in Eq. ([3]) are truncated, due to the lack 
of exact Hermiticity when such approximations are made. However, it is clear that the 
important Hellmann-Feynman theorem is preserved in all such approximations. 
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We also note that Eq. (^) represents a coupled set of nonlinear multinomial equations 
for the c- number correlation coefficients {Si}. The nested commutator expansion of the 
similarity-transformed Hamiltonian, 

H = e~ s He s = H+[H,S] + ±[[H, S], S] + ■ ■ ■ , (11) 

together with the fact that all of the individual components of S in the sum in Eq. 
commute with one another, imply that each element of S in Eq. (§) is linked directly to 
the Hamiltonian in each of the terms in Eq. (0). Thus, each of the coupled equations 
(|) is of linked cluster type. Furthermore, each of these equations is of finite length when 
expanded, since the otherwise infinite series of Eq. ([□]) will always terminate at a finite 
order, provided (as is usually the case) that each term in the second-quantised form of the 
Hamiltonian H contains a finite number of single-body destruction operators, defined with 
respect to the reference (vacuum) state |$). Therefore, the CCM parametrisation naturally 
leads to a workable scheme which can be efficiently implemented computationally. It is also 
important to note that at the heart of the CCM lies a similarity transformation, in contrast 
with the unitary transformation in a standard variational formulation in which the bra state 
(\&| is simply taken as the explicit Hermitian adjoint of \^). 

The CCM formalism is exact in the limit of inclusion of all possible multi-spin cluster 
correlations for S and S, although in any real application this is usually impossible to 
achieve. It is therefore necessary to utilise various approximation schemes within S and 
S. The three most commonly employed schemes previously utilised have been: (1) the 
SUBn scheme, in which all correlations involving only n or fewer spins are retained, but 
no further restriction is made concerning their spatial separation on the lattice; (2) the 
SUBn-m sub- approximation, in which all SUBn correlations spanning a range of no more 
than m adjacent lattice sites are retained; and (3) the localised LSUBm scheme, in which all 
multi-spin correlations over distinct locales on the lattice defined by m or fewer contiguous 
sites are retained. The specific application of the CCM to the spin-half transverse Ising 
model in the paramagnetic and ferromagnetic regimes is now described. 
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1. The Paramagnetic Regime 



In the paramagnetic regime, a model state is utilised in which all spins point along the 
x-axis, although it is found to be useful to rotate the local spin coordinates of these spins 
such that all spins in the model state point in the 'downwards' direction (i.e., along the 
negative z-axis). This (canonical) transformation is given by, 

a x -> -a z , a v ^ a y , a z ^ a x , (12) 

such that the transverse Ising Hamiltonian of Eq. ([!]) is now given in the (rotated) spin- 
coordinate frame by, 

h=(~ + \)n-y, Um + *r°j + °7°t + °t°j\ + a £ °* > (!3) 

VZ 7 (id) 1 J i 

where = \{a% ± icr^) . In these local coordinates the model state is thus the "ferromag- 
netic" state |^) = | 11 ■ ■ ■ I ■ ■ ■) in which all spins point in the downwards direction. In 
order to reflect the symmetries of this Hamiltonian, the cluster correlations within S are 
explicitly restricted to those for which = J2i sf (in the rotated coordinate frame) is an 
even number. Hence, the LSUB2 approximation is defined by 

^|l^> ( 14 ) 

where p covers all nearest-neighbour lattice vectors. The ground-state energy is now given 
in terms of b\ by, 

— = -(1 - 60 . (15) 

It is found that this expression is valid for any level of approximation for S. Using Eq. (0) 
it is found that, 

5bl + 4A&! -1=0, (16) 

and hence an approximate solution for the ground-state energy at the LSUB2 approximation 
level purely in terms of A may be obtained. The SUB2 approximation contains all possible 
two-body correlations, for a given lattice, and is defined by 
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where the index r indicates a lattice vector. Eq. (0) may once again be utilised to determine 
the SUB2 ket-state equations. Hence the CCM SUB2 ket-state equation corresponding to a 
two-body correlation characterised by index s is given by, 

(1 + 2bl)5 StP - 46A + £ b r b r+s+p + 2b s+p - AXb s = . (18) 

r 

We note that this equation is meaningful only for s ^ as we may only ever have one 
Pauli raising operator per lattice site. This equation may be solved by performing a Fourier 
transformation. (Details of how this is achieved in practice are not given here and the 
interested reader is referred to Refs. 0, p3|.) An alternative approach, however, is to use 
Eq. ( |T8| ) in order to fully define the SUB2-m equations. This is achieved by truncating 
the range of the two-body correlations (i.e., by setting \s\ < m), and the corresponding 
SUB2-m equations may be solved numerically via the Newton-Raphson technique (or other 
such techniques). We note that coupled sets of high-order LSUBm equations may be derived 



using computer-algebraic techniques, as discussed in Ref. [24]. The technicalities of these 



calculations are not considered here, but the interested reader is referred to Ref. |p4j| . A full 
discussion of the CCM results based on the paramagnetic model state is deferred until Sec. 
III. 



2. The Ferromagnetic Regime 

In the ferromagnetic regime, a model state is chosen in which all spins point 'downwards' 
(along the negative z-axis), and so the Hamiltonian of Eq. ([I]) may therefore be utilised 
directly within the CCM calculations. The lowest order approximation is the now SUB1 
approximation (in which case, S = aj^i&t) an d the ground-state energy is given in terms 
of a by, 

f = A(l-o) • (19) 
It is again noted that this expression is valid for any level of approximation in S. In this 
case, it is found that the solution of the SUB2 approximation collapses onto the LSUB2 
solution due to the simple nature of the Hamiltonian and model state, although it is again 
possible to perform high-order LSUBm calculations. Furthermore, the lattice magnetisation 
(i.e., the magnetisation in the z-direction) , M, is defined within the CCM framework by, 

1 N ~ 

M = | , (20) 

i=l 



which may be determined once both the ket- and bra-state equations have been solved at 
a given level of approximation. Again, the discussion of the results for this model state is 
deferred until Section III. 



B. The CBF Formalism 

The treatment of the transverse Ising model by the CBF method is begun by defining 
the lattice magnetisation (i.e., again the magnetisation in the z-direction), given by 

" = (21) 

for a ground-state trial wave function, \vp). Furthermore, the 'transverse' magnetisation is 
given by, 

_ (y> | of | j,) 

(m ■ (22) 

It is also found to be useful to define a spatial distribution function (which plays a crucial 
part in any CBF calculation) in the following manner, 

(V> | <7f<7 2 | V) 

#( n ) = 7Tii\ ' ( 23 ) 



where n = — Tj. The corresponding approximation to the ground-state energy per spin is 
given by 

#-^-|- ? A(^ + A(l-^. (24) 

where the function A(n) is equal to unity when n is a nearest-neighbour lattice vector 
and is zero elsewhere. It is noted that the distribution function g(n) may be decomposed 
according to g(n) = 5 nfi + (1 — 5 n fi)M 2 + (1 — M 2 )G(n) such that G(n) now contains the 
short-range part of the spatial distribution function and vanishes in the limit |n| — > oo. 
The magnetisation, M, and the transverse magnetisation, A, may now be expressed in a 
factorised form in terms of a 'spin-exchange strength', riu, such that, 

A = (1 - M 2 )^n 12 . (25) 

The energy functional is now expressed in terms of G(n) and n i2 as, 

^ = (l-^-EA(n)G(n)} + A{l-(l-M 2 )%} . (26) 
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Note that in the mean-field approximation G(n) in Eq. ( p6[) is set to zero (for all n) and 
rii2 is set to unity. 

In order to determine the ground-state energy and other such ground-state expectation 
values, a Hartree-Jastrow Ansatz is now introduced, given by 

\^}=exp{MU M + U}\0} . (27) 

The reference state |0) is a tensor product of spin states which have eigenvalues of +1 with 
respect to a x . The correlation operators U and Um are written in terms of pseudopotentials, 
n(r.y), 1*1(1-4), and u M (r i:j ), where 

1 N 

i<j 

and 

N ^ N 

U M = E«i( r *K + ^EM^IW + of) . (29) 

The pseudopotential «i(rj) = m x is independent of the lattice position by translational in- 
variance, and the pseudopotentials, w(ry) and UM(jij), similarly depend only on the relative 
distance, |n| = | — | = |r^ |. The Jastrow correlations are determined via a cluster ex- 



pansion of the various quantities in the Hamiltonian, as explained in Refs. p9| , p0| , pT| , |32 
A common approximation yields the hypernetted chain (HNC/0) equations, which one may 
solve iteratively in order to determine the Hartree-Jastrow pseudopotentials. One then 
wishes to determine the expectation values such as the ground-state energy, and in the 
paramagnetic regime an explicit assumption is made that M = 0. However, in the ferro- 
magnetic regime M is taken to be a variational parameter with respect to the ground-state 
energy of Eq. (^6|) . 

There are now two ways of determining the pseudopotentials from the HNC equations. 
The first such approach is to assume that the pseudopotential has the simple parametrised 
form, 

u(n) = a A(n) , (30) 

where A(n) is unity if n is a nearest-neighbour lattice vector and is zero otherwise. This 
approach is henceforth denoted as the parametrised HNC CBF method. The c-number 
a is taken to be a variational parameter with respect to which the ground-state energy is 
minimized. In the paramagnetic regime, the minimum of the energy surface as a function of a 



is sought, at a given value of A. This is easily performed computationally, and the solution is 
readily tracked iteratively, starting from the trivial limit A — ► oo and then moving to smaller 
values of A. In the ferromagnetic regime, one again searches for a minimum of the energy 
surface, but this time with respect to both a and M, at a given value of A. In this case, one 
tracks from the trivial limit of A = to higher values of A. In previous articles]^, [JI], [JT], |3"2"| ], 
this was achieved by analytically determining the derivative of the energy with respect to M, 
although in this article a computational minimisation of the energy is performed with respect 
to both variables. The second such approach allows one to find the optimal pseudopotential 
within the CBF/HNC framework from a functional minimisation, 



Within the context of this article, this is henceforth denoted as the paired phonon approxima- 
tion (PPA) (or, more precisely, the paired magnon approximation), and the corresponding 
equations are the PPA equations. Note that in this article the PPA calculation is only 
performed in the paramagnetic regime, although PPA results in the ferromagnetic regime 
have also been performed previously p0|, plf| . (Indeed, results for the phase transition points 
predicted by the PPA CBF approach of Ref.JSy, [JT| are quoted in the Table |.) A full 
discussion of the results of the CBF calculations presented in this section, in comparison to 
the corresponding results of the CCM and the VQMC method, is given in Section III. 

C. The VQMC Formalism 

Although the specific variational calculations presented in this article concentrate on the 
spin-half transverse Ising model, we note that the formalism presented in this section is given 
in a generalised form and that the treatment of other spin models would follow a similar 
pattern. We shall specifically consider here the transverse Ising model in the ferromagnetic 
regime where the relevant Hamiltonian is defined by Eq. ([[]). An Ansatz for the expansion 
coefficients, {c/}, of a ground-state wave function defined by 



is chosen. Note that {|/)} denotes a complete set of Ising basis states, defined as all possible 
tensor products of states on all sites having eigenvalues ±1 with respect to o z . An expression 




(32) 



i 
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for the ground-state energy is thus given by, 

E = Y, h ,i 2 c* h c l2 {h\H\I 2 ) ^ 

J2r \ c i'\ 2 

Specifically for the spin-half transverse Ising model, a Hartree-Jastrow Ansatz p0| (for 
A > 0) is now defined with respect to the {c/} expansion coefficients, where 

Cl =(I\l[(l + a k Pl)]l(l + fJpJpj + Plpj})\I) . (34) 
k v y i<j v L J/ 

The P^ and P^ are the usual projection operators of the spin-half 'up' and 'down' states 
respectively. The simplest form of the variational Ansatz of Eq. (133) is given by, 



fi if % and j are nearest neighbours, 
/y= { (35) 
otherwise. 

The symmetry-breaking term a^(= a) is also independent of i by translational invariance. 
The expectation value of Eq. (^) may now be evaluated directly and the variational ground- 
state energy minimised with respect to both a and fi at each value of A. However, such a 
calculation is soon limited by the rapidly increasing set of Ising states and the amount of 
computational power available. Indeed, for the spin-half transverse Ising model the number 
of states that one must sum over is 2^, where N is the number of lattice sites. For the linear 
chain it is possible to solve for chains of length N ~ 12 with relatively little computational 
difficulty, although the calculations with N > 12 grow rapidly in computational difficulty 

Hence, as an alternative for lattices of larger size, we may simulate the summation over 
all the states I\ in Eq. fl3"3"|). In order to to do this we define the probability distribution for 
the set of states {|/)}, given by 

P(I) - ^~r 2 , (36) 
and the local energy of these states, given by 



W=E-(Wi) • (37) 
h Cl 



The expression of Eq. (|3~3D may thus be equivalently written as, 



E = J2 P{I) E L (I) . (38) 
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We now wish to perform a random walk based on the probability distribution of Eq. (j36[) . 
However, a few more useful quantities are best defined before a detailed description of the 
VQMC algorithm is actually given. Firstly, the acceptance probability, A(I — > /'), of Monte 
Carlo 'move' from state \I) to state \I') is given by 

A(I -> J') = min [1, g(J -> /')] , (39) 

where 

g (j _> /') = P ( J/ ) T ( J/ ~> J ) . (40) 

yi J P(/)T(J -> J') 1 J 

Here, T(J — > /') is the sampling distribution function. For spin lattice problems, if state 1 1) 
can connect to K(I) possible Ising states via the off-diagonal elements of the Hamiltonian 
then T(I — > I') = 1/K(I). Hence, q(I — > I') is written as 

P{I')K{I) 

qiI ^ I) = P(I)K(I>) • (41) 
For the spin-half transverse Ising model, we note that K(I) is equal to iV for any state |J) 
and so the common factor of N in Eq. (|4l| ) cancels. The simplest VQMC procedure is now 
defined by the following algorithm: 

1. Select an initial Ising state \I) for which ^ 0, where |^) is the 'true' ground-state 
wave function of the system. 

2. Choose a particular state |J') out of the K(I) possible states accessible to |J) via the 
off-diagonal elements in H. 

3. Define a random number r in the range [0, 1] and accept this move from state \I) to 
state |/') if and only if 

A(I ->/')> r . (42) 

4. If the move is accepted then let I — > I' and cj — > cp. 

5. Obtain the local energy E(I) of Eq. (|37j) for state |J). 

6. Repeat from stage (2) N MC times. 

7. The average ground-state energy (and the error therein) may be determined from the 
Nmc number of local energies during the simulation. 
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The minimal VQMC ground-state energy is now obtained by searching over the variational 
parameter space for either the lowest ground-state energy or lowest variance in the ground- 
state energy (here for a given value of A in if). In order to determine the lattice magneti- 
sation, we note that 



M = 



i 

where P(I) is the probability density function given above, and 



h 



N 



(43) 



(44) 



Hence, a mean value (and its associated error) for the VQMC lattice magnetisation, M, 
may be obtained by determining the average of the local lattice magnetisation, M L (I), 
throughout the lifetime of the run. A discussion of the results of the variational calculations 
discussed here is given in Section III. 



D. The Infinite Lattice Limit and Convergence of Results 

In this section we consider how the results of each method are determined in the infinite 
lattice limit. Firstly, it is noted that the CCM method produces expectation values which 
are size-extensive (i.e., the numerical values of each expectation value scale linearly with 
N), and we always deal with an infinite lattice in all calculations from the very outset. 
Furthermore, the 'raw' CCM LSUBm results based on the ferromagnetic model state are 
found to converge rapidly with increasing LSUBm approximation level m over most of the 
ferromagnetic regime except for a region very near to the phase transition point. In order 
to obtain even better results for the CCM method across the whole of this regime, a simple 
extrapolation of the LSUBm data in the limit m — > oo has also been carried out at each 
value of A separately. The results of the extrapolation using a leading-order 'power-law' 
dependence (see Appendix A for details) are denoted as Extrapolated(l) CCM results and 
results of a Pade approximation for / = (again, see Appendix A for details) are called 
Extrapolated(2) CCM results. In the paramagnetic regime, the results for the ground-state 
energy are found to converge extremely rapidly with LSUBm approximation level over the 
whole of this regime and so no extrapolation of these results is necessary. 
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For the CBF method, although the treatment presented here is formally valid for any 
lattice size (including the infinite lattice case), the results presented have been obtained for 
finite-sized lattices. The results are found to converge extremely rapidly with increasing 
lattice size, and the results of the 20x20 square lattice (used in the figures given below) are 
found to be essentially fully converged for all A except for a region very near to the critical 
point. 

The results for the VQMC method presented below have been obtained for a 16x16 
square lattice, where the number of Monte Carlo iterations was set to 50000. As for the 
CBF calculations, the 16x16 lattice is again expected to be large enough for the VQMC 
results to be essentially fully converged to the infinite lattice limit for all values of A except 
for a region near to the critical point. By comparing the results of a 10x10 lattice VQMC 
calculation with those of a 16 x 16 lattice it was found that this was indeed true. Furthermore, 
we note that for small A the variational minimum of the ground-state energy was found to be 
rather flat with respect to the variational parameters a and fx, and the ground-state energy 
was also highly converged with increasing lattice size. However, as the phase transition point 
is approached one finds that a precise evaluation of the position of the variational minimum 
with respect to a and fx becomes harder to determine. 

III. RESULTS 

The results for the ground-state energy per spin of the spin-half transverse Ising model 
on the square lattice in the ferromagnetic regime are shown in Fig. [I]. We see from this 
figure that excellent correspondence between the results of the CCM, CBF, and VQMC 
methods for the ground-state energy is obtained in this regime. We note, however, that the 
extrapolated CCM results appear to lie very slightly lower than the other two sets of results, 
especially near to the phase transition point. This indicates the increasing importance of 
higher-order correlations for the ground-state energy near to the phase transition point. 

The results for the lattice magnetisation M obtained using the CCM, CBF, and VQMC 
formalisms are shown in Fig. ^| (and also Fig. |] in Appendix A) for the spin-half transverse 
Ising model on the square lattice. We note that the 'raw' CCM LSUBm results for the lattice 
magnetisation do not become zero at any value of A for any finite value of the truncation 
index m, because all of the ferromagnetic order inherent in the model state must be destroyed 



14 




FIG. 1: Results for the ground-state energy per spin, E g /N, of the spin-half transverse Ising model 
on the square lattice in the ferromagnetic regime using the CCM, CBF, and VQMC approaches. 

in order for M to be zero in this case. In practice this is a difficult thing for the CCM to 
achieve with this model state. However, we may see from Fig. ^| that the extrapolated 
CCM results are in excellent agreement with those results of the CBF and VQMC methods. 
In addition, it is possible to imagine other CCM model states, such as a spin-flop model 
state or a mean-field model state, in which the lattice magnetisation with respect to this 
state is not a priori fully saturated for all A. Furthermore, we remark that a treatment of 
this problem with the ferromagnetic model state using the extended coupled cluster method 
(ECCM) [|], [5], [7], ^8] would present an interesting challenge. 



In the paramagnetic regime, the ground-state energy per spin for the transverse Ising 
model on the square lattice is presented in Fig. |3|. The 'raw' CCM LSUBm results for the 
ground-state energy are already highly converged with increasing truncation index m, even 
up to the phase transition point. We note again that an extrapolation in the limit m — ► oo 
is therefore not necessary. Indeed, good correspondence between the results of the different 
methods is seen although it is noted that CCM LSUB4 and LSUB6 ground-state energies 
lie lower than those predicted by the CBF. This indicates that high-order order correlations 
become increasingly important the nearer one gets to the phase transition point. However, 
this could be rectified, in principle, for the CBF method by the inclusion of higher-order 
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FIG. 2: Results for the lattice magnetisation, M, of the spin-half transverse Ising model on the 
square lattice using the CCM, CBF, and VQMC approaches. 

(than pairwise Jastrow) correlations in the ground-state wave function. We note that in 
practice, however, the inclusion of such higher-order correlations in the ground-state wave 
function in the CBF and VQMC methods is a difficult and unresolved question. 

It is also possible to determine the second-derivative of the ground-state energy per spin 
with respect to A for the CCM calculations based on the paramagnetic model state, defined 



It is found that x diverges at some critical value A c for the SUB2 approximation in any 
dimension, and for SUB2-m and LSUBm (with m > 4) approximations for spatial dimen- 
sionality greater than one. Again, this behaviour is associated with a phase transition in 
the real system and the point at which this occurs is denoted A c . Correspondingly, CCM 
results for A < A c based on the paramagnetic model state do no exist, and hence A c acts as a 
terminating point for the calculation in the paramagnetic regime. Also, it is found that the 
SUB2-m results for A c as a function of m scale with m~ 2 and a simple linear extrapolation 
gives the full SUB2 result for the critical point to within a 2% accuracy. By analogy, this 
rule has also been used for the LSUBm results to extrapolate to the limit m — > oo, and the 
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FIG. 3: Results for the ground-state energy per spin, E g /N, of the spin-half transverse Ising model 
on the square lattice in the paramagnetic regime using the CCM and CBF approaches. The boxes 
for the CCM data indicate the terminating points at which \ becomes infinite. 

results thus determined are shown in Table |. The values thus obtained are also in good 
agreement with the points at which M — > from the extrapolations in the ferromagnetic 
regime discussed above (and see Fig. 2 for the square lattice case). 

For the CBF and VQMC methods the point, in terms of A, at which M becomes zero 
is taken to indicate a quantum phase transition and is again denoted, A c , and these results 
are presented in Table |. The phase transition point predicted by the VQMC method on 
the square lattice case is estimated to be at A c = 3.15 ± 0.05. For the linear chain, the 
expression for ground-state energy of Eq. ( p3|) has been obtained directly for chains with 
N < 12. These results are found to be in good agreement with a previous calculation pl| 
using the Ansatz of Eq. ([34]) for the linear chain transverse Ising model which predicted a 
value for the phase transition point of A c = 1.206. 
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TABLE I: Results for the critical points of the spin-half transverse Ising model on various lattices 
using the CCM and CBF approaches. These results are compared to those of exact calculations for 
the linear chain [46] and to RPA |^7| and exact cumulant series expansions [49| in higher spatial 
dimensionality. For the linear chain, we note that the LSUBm approximation does not show any 
evidence of a critical point, at the levels of approximation shown, and this is indicated by 'none'. 
(Previous CBF results of Refs. |2^, |30|, [3l]] for the square and cubic lattices are also appropriately 
indicated.) 





Linear Chain 
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Triangular 


Cubic 


Classical 
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RPA a 




3.66 




5.76 


CCM A c SUB2 


1.44 


3.51 


5.42 


5.55 


CCM A c LSUB4 


none 


2.41 


4.27 


3.85 


CCM A c LSUB6 


none 


2.76 


4.57 


4.61 


CCM A c LSUBcx) 


none 


3.04 


4.81 


5.22 


parametrised HNC CBF 


1.22 


3.12 6 


4.91 


5.17(6/ 


PPA CBF C 




3.14 




5.10 


Variational or VQMC Calculations 


1.206 d 


3.15±0.05 






Exact or Series Expansions Calculations 


1.0 e 


3.044/ 


4.768[/y 


5.1531// 



°from Ref. |4J 
6 from Ref. @ 
c from Refs. §(J |l 
d from Ref. @ 
e from Ref. @ 
/from Ref. M 



IV. CONCLUSIONS 



In this article, results of the CCM, CBF, and VQMC approaches for the ground-state 
energy, the lattice magnetisation, and the position of the critical point of the spin-half 
transverse Ising model on various lattices have been presented. These results have been 
seen to be in excellent agreement both with each other and with those results of exact 
calculations for the linear chain and those of exact cumulant series expansions for higher 
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spatial dimensionality. Indeed, by treating these systems using three separate approaches, 
it has been shown that each set of results has been mutually supported and reinforced by 
those of the other approaches. 

Furthermore, we have gained some insight into the strengths and weaknesses of each 
approach. This is exemplified in the different parametrisations of the ground-state wave 
function. The CBF and VQMC approaches both utilise Jastrow wave functions and their 
bra states are always the explicit Hermitian adjoints of the corresponding ket states. Hence, 
for the CBF and VQMC approaches, an upper bound to the true ground-state energy is, 
in principle, obtainable, although the approximations made in calculating the energy may 
destroy it. By contrast, the CCM uses a bi- variational approach in which the bra and 
ket states are not manifestly constrained to be Hermitian adjoints and hence an upper 
bound to the true ground-state energy is not necessarily obtained. Also, the CCM uses 
creation operators with respect to some suitably normalised model state in order to span 
the complete set of (here) Ising states. The other approaches, in essence, use projection 
operators to form the Hartree and the Jastrow correlations. For the CBF case, this is with 
respect to a reference state, whereas for the VQMC case, the Hartree- Jastrow Ansatz is 
encoded within the expansion coefficients of the ground-state wave function with respect to 
a complete set of Ising states. In some sense, the CCM is found to contain less correlations 
than the others at 'equivalent' levels of approximation (e.g., the CCM LSUB2 approximation 
versus Hartree and nearest-neighbour Jastrow correlations). A fuller account of the different 
parametrisations of the ground-state wave function within the CCM and CBF methods has 
been given in Ref. ||32[| . However, in practice the other methods are difficult to extend to 
approximations which contain more than two-body or three-body correlations. By contrast, 
the CCM is well-suited to treat such higher-order correlations via computational techniques, 
as has been demonstrated here. Furthermore, the CCM requires no information other than 
the approximation in S and S in order to determine an approximate ground state of a given 
system. The CBF method, however, may require that only a certain subset of all possible 
diagrams are summed over (e.g., the HNC/0 approximation). The VQMC approach also 
often requires an intimate knowledge of the manner in which the two-body correlations 
behave with increasing lattice separation if all two-body correlations are to be included. 
This information may be approximated, for example, by use of the results of spin-wave 
theory. In any case, it is often necessary to reduce the minimisation of the variational 
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ground-state energy with respect to iV parameters to much fewer parameters. Another 
potential application of all of the methods presented here is the use of their ground-state 
wave functions as trial or guiding wave functions in (Green function or similar) quantum 
Monte Carlo calculations. 

Finally, encouraged by these results for the transverse Ising model, we intend to extend 
them to other models of interest, such as systems with higher quantum spin number or those 
with complex crystallographic lattices. A further goal is to extend the treatment of this and 
other spin models, via these methods, to non-zero temperatures. 
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APPENDIX A: EXTRAPOLATION OF CCM RESULTS 



In this Appendix, we explain how we extrapolate a set of LSUBm data points, {xi,yi}, 
in the limit i — > oo at each value of some parameter (A) within the Hamiltonian separately. 
Note that the number of data elements to be extrapolated is given by the index p. The 
value of Xi is now set to be 1/m and yi is set to be the corresponding value of an expectation 
value (for example, the lattice magnetisation) determined using the CCM at this level of 
approximation at a given value of A. Note that the value of m must increase with increasing 
index i, although m and % do not have to be equal. 

Before the extrapolation procedures are given in detail, we define some useful quantities. 
Firstly, the mean value of a set {q} is denoted by c and of a set {di} is denoted by d. 
Secondly, the linear correlation, R, of a set of two-dimensional points, {q, di}, is defined by 

r = - £?=ite- g )(*-3) . (A1) 
'TLiici-cFJTLiidi-d)* 



We are now in a position to outline the the first extrapolation procedure. This procedure 
firstly assumes that the data scales with a leading-order "power-law" dependence, given by 

y i = a + bx v i . (A2) 

We set Q = log(xj) and di = \og{yi — a), where {xi, yi\ is the LSUBm data set at some fixed 
value of a parameter within the Hamiltonian. Hence the best fit of the data set, {xi,yi}, to 
the power-law dependence of Eq. ( |A2| ) is obtained when the absolute value of R is maximised 
with respect to the variable a. Indeed, we make the assumption that this value of a is then 
taken to be the extrapolated value of the in the limit i — > oo (in which case, Xi — > 0). 

The second extrapolation procedure of the LSUBm data uses Pade approximants. This 
is achieved by firstly assuming that the set of data can be modeled by the ratio of two 
polynomials, given by 

1 + E j= i bjxi 

Note that when I = 0, this is a simple integral power series. This furthermore implies that, 

a + a x Xi + a 2 x1 H h a k x\ = y i J r (hxi + b 2 x\ H V bix\)yi . (A4) 

We now wish to determine the coefficients a,- and bj in order to find the polynomials in Eq. 
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, and Eq. 



is rewritten in terms of a matrix given by, 



^ 1, X\, xf, 
1, X2, X21 



•, x 1 

k 



\ 



-61 



2/2 



(A5) 



The inverse of the matrix in Eq. 
determined. (Note that k + l + 1 = p.) However, we also note that because Xi — > as i — > 00, 
ao gives us the extrapolated value of {t/i} in the limit m — > 00. Furthermore, using this 



is now obtained and the coefficients dj and bj are 



method with I = we found that a previous extrapolated result [^4] of CCM LSUBm data 
for the sublattice magnetisation of the Heisenberg antiferromagnet on the square lattice was 



reproduced. In this previous calculation the sublattice magnetisation was extrapolated 
in the limit m —>■ 00 by fitting the LSUBm points (with m = 4, 6, 8) to a quadratic function 
in 1/m thus giving an extrapolated value of about 0.6. 

In this article, we have already plotted extrapolated CCM LSUBm results for the lattice 
magnetisation of the square lattice spin-half transverse Ising model in Fig. [| and these 
results were seen to be in excellent agreement with those results of CBF and VQMC calcu- 
lations. However, a further discussion of the extrapolated CCM results presented here is also 
useful in order to illustrate the strengths and weaknesses of the extrapolation procedures 
outlined in this appendix. We can see from Fig. [| below that the results for the Pade ap- 
proximant extrapolation with I = 3 contains a zero in the denominator of Eq. (|A3|) at about 
A ~ 2.6 such that the results show a divergence for M in Fig. |] which is simply an artifact 
of the extrapolation procedure. This is because an assumption is made as to the scaling 
of the LSUBm data with 1/m to some functional form. The validity of this assumption is 
unknown as no exact scaling laws are known, as yet, for the behaviour of CCM LSUBm 
results as functions of m. However, the empirical evidence in Fig. [| suggests that this is 
a reasonable assumption over much of the ferromagnetic phase, except, of course, for those 
points at which the Pade approximant results demonstrate this 'artificial' divergence. 
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FIG. 4: Results for the lattice magnetisation, M, of the spin-half transverse Ising model on 
the square lattice in the ferromagnetic regime for the CCM LSUBm approximation with m = 
{2, 3,4, 5, 6, 7}. The extrapolation of the LSUBm results at each separate value of A is performed 
in two ways. The first uses Pade approximants in order to perform the extrapolation, and the 
second assumes a leading-order power-law scaling of the lattice magnetisation with vrT 1 . 
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